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Abstract. Imaging of perfectly conducting crack(s) in a 2-D homogeneous medium using 
boundary data is studied. Based on the singular structure of the Multi-Static Response (MSR) 
matrix whose elements are normalized by an adequate test function at several frequencies, an 
imaging functional is introduced and analyzed. A non-iterative imaging procedure is proposed. 
Numerical experiments from noisy synthetic data show that acceptable images of single and 
multiple cracks are obtained. 



1. Introduction 

One considers the inverse scattering problem of the retrieval of the shape of perfectly conducting 
cracks buried within a homogeneous medium. Data are scattered field measured at the search 
domain boundary. 

To our knowledge, inverse scattering with a single, perfectly conducting crack in 2-D satisfying 
a Dirichlet boundary condition has been studied in [3]. In this reference, the retrieval algorithm 
was based on a Newton-type iterative scheme so that a good initial guess close to the unknown 
crack was needed to ensure convergence. But an iterative scheme generally requires a large 
amount of computational time, and it is hard to extend to the reconstruction of multiple cracks 
unless strong hypotheses are made beforehand. 

So, non- iterative reconstruction algorithms have been developed, e.g., Multiple Signal 
Classification (MUSIC)-type ones [21 El ED] 5 topological-derivative-based ones [HIE]. Recently, a 
multi-frequency algorithm has been proposed; related works are in [H [61 [11] and references. 
But in many cases, the analysis is not fully reliable, in particular since phenomena such as the 
unexpected appearance of ghost replicas cannot be explained. 

Here, one is proposing an improved imaging algorithm for perfectly conducting cracks from 
scattered field data measured at the search domain boundary, and one analyzes its structure 
by designing a multi-frequency imaging functional weighted by each frequency. This involves 
that elements of the Multi-Static Response (MSR) matrix normalized by a test function can 
be represented in a suitable asymptotic form. Via careful analysis, one can infer a relationship 
between the proposed imaging functional and Bessel functions of integer order and explain why 
this imaging functional yields location and shape of the cracks. 

The discussion goes as follows. In section [21 the direct scattering problem is summarized and 
the asymptotic expansion formula introduced. In section [3] the multi-frequency-based imaging 



functional weighted by each frequency is put forth and its properties discussed. Numerical 
experiments are in section 0] for single and multiple cracks. Conclusion is in section [5j 



2. Direct scattering problem and asymptotic expansion formula 

In this section, one considers 2-D electromagnetic scattering by a perfectly conducting crack, 
r, fully buried within a homogeneous domain Si C K 2 . One assumes that T is a narrow, linear 
crack represented as follows 

r = {z : z = (x,y),-p < x < p} , (1) 

and that it does not touch the boundary <9Sl, i.e., T n <9Sl = 0. 

Throughout this contribution, we remain with the Transverse Magnetic (TE) polarization 
case: at a given frequency u>k ^> p, let u( ra )(x;c<jfc) be the time-harmonic total field that satisfies 

' AuW (x ; u k ) + oj 2 k u^ (x; uj k ) = for x G fi\F, 

n( n )(f;w fe )=0 for x G T, 

9u(")(s;a; fc ) ae^ 



dv(x) dv{x) 



for x G 90. 



Here, z? is the outward unit normal vector to <9Sl and {6 n : n = 1, 2, • • • , iV} are the iV vectors 
that describe the directions of incidence, assumed equi-distributed on the unit circle § . In the 
same manner, ujj (x;u;&) satisfies ([2]) without T. Then, the following asymptotic relationship of 
u^(x;0Jk) and Uq (x;wfe) holds uniformly on x G <9Sl (see [2] for instance): 

n (n) (x;w fc ) - u[, n) (x;u; fc ) = m ^ 2 ) ^ & z;u k ) + of - j , (3) 

where < p < 2 and A/"(x, z; uj k ) is the Neumann function for Si, 

AA/"(x, z; Wfe) + u)fcN(x, z;cok) = —8(x, z) for fed, 

dN%z;u k ) (4) 

= for x G ai 2. 

a^(x) 

3. Non-iterative multi-frequency imaging algorithm from MSR matrix 

One assumes that there exist L different small cracks Ti with same length p centered at z*i, 
I = 1, 2, • • • , L, in Si. We apply the asymptotic formula ([3]) to establish an imaging functional of 

the crack. For this, one uses the eigenvalue structure of the MSR matrix M(u;fc) = yMnm 

(k) 

where Mnm denotes the normalized boundary measurement data: 



(fc)^ 



n,m=l 



Here, one eliminates the residue term of ([3]) and chooses a test function v^ m \x;ujk) = e~ lu)kdm ' x . 

Now, one designs an imaging functional. First, one performs the Singular Value 
Decomposition (SVD) of M(wfc), so that 



— T 

M(w fc ) « J3«7iK) ff iK)^i H), 



where &i(uj k ) are (non-zero) singular values and where Ui(co k ) and Vi(oj k ) are left and right 
singular vectors. Then, one introduces a weighted imaging functional at several discrete 
frequencies {u) k : k = 1, 2, • • • , K} as 

E(x;K) := £ • £/,H)) w fc ) ■ ?,K)), (6) 

z — i ;— i \ / \ / 
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where 
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W(x; u k ) = " ' " with £>(£; Wfe ) = e l ^ x , e 
\D(x;u k )\ V 

It should be noted that W(x;u k ) satisfies the following properties (see [I]) 

Ul(u) k ) ~ W(zi;uj k ) and ~ W^J^*.)- (7) 

If one assumes that numbers N and K are large enough, applying (JJJ) to dU) yields 
L K / N \ / N \ 

E(x; if ) ~ £ I>M E e*"*^- 50 ) ( X ) 
1=1 fc=l ^ s=i / V t=i ' 

1 L 

i=l 



^ (joO^n) 2 + Ja(o;xn) 2 ) " y (^in) 2 + Ji^nA 



where J p (-) denotes the Bessel function of order p and of the first kind. A full derivation is to 
appear in an extended version of this contribution. Some properties of E(x; K) are as follows: 

• Jq{oj\zi — x\) 2 + J\{u)\zi — x\) 2 reaches its maximum value at oj\zi — x\ = 0. Hence, |E(x; K)\ 
plots are of large and small magnitudes at x = z\ G T\ and at x ^ ^ £ Tj, respectively. 

• When K is small, |E(x; K)\ plots exhibit replicas at local maxima and minima of Jo- 

• ./V must be large enough. It means that if one applies this algorithm with a small N, 
poor results follow. Moreover, by looking at ([5]), this analysis cannot be applied to the 
limited-view inverse scattering problem. 

• Using frequencies high enough guarantees successful imaging performance. 

• Based on the Statistical Hypothesis Testing pQ, the proposed algorithm is robust vs. random 
noise when K is large enough. 

4. Numerical experiments 

In this section, one introduces numerical examples for imaging small and long arc-like cracks. 
The operation frequency is taken of the form uo k = 2ir/X k ; here A^, k = 1,2, ■■■ , if(= 10), 
is the given wavelength. The uo k are equi-distributed in the interval [g^wk-] with Ai = 0.6 
and \k = 0.4. A white Gaussian noise with 20dB signal-to-noise ratio (SNR) is added to the 
unperturbed data. 

Figure [1] displays maps of |E(x; 10) | for three small cracks (left), a long single crack (center), 
and two long cracks close to one another (right). Notice that the analysis above is restricted 
to the small crack case, but it can be applied to extended (long) ones. Although some replicas 
appear since the chosen value of K is rather small (= 10), |E(s?;if)| plots exhibit maximum 
amplitudes at the location of the cracks so that one can recognize their shape. Taking them 
as initial guess could yield more accurate shapes via a Newton-type algorithm [4J or a level-set 
method [3j. 
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(a) \E(x; 10) | with N = 12 




(b) \E(x; 10)| with N = 32 



(c) |E(x; 10)| with N = 36 



Figure 1. Map of |E(»;if)| with K = 10 for various cracks. 



5. Conclusion 

Based on the structure of the singular vectors of the MSR matrix at several discrete frequencies of 
operation, an improved image function motivated from Kirchhoff migration has been proposed 
to image perfectly conducting cracks. The proposed imaging functional is related to Bessel 
functions of integer order. Numerical experiments with noisy data are discussed and they 
show the effectiveness of the corresponding imaging algorithm. The current contribution deals 
with a full-view inverse scattering problem. Application to a limited-view inverse problem is 
forthcoming work. Along the same line of thought, extension to cracks with Neumann boundary 
condition should be an interesting subject. 
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